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INVERSE  SCATTERING  THEORY  AND 
ALMOST  PERIODIC  MEDIA 


I.  INTRODUCTION 

A  fundamental  problem  in  physics  is  the  investigation  and  description  of  an  unknown  object  by 
probing  or  sensing  with  electromagnetic  waves.  Inverse  scattering  theory  is  concerned  with  the  analyti¬ 
cal  and  computational  methods  for  solving  this  general  problem.  The  physical  model  that  will  be  con¬ 
sidered  in  this  report  is  the  determination  of  the  effective  permittivity  and  density  distribution  of  an 
inhoniogeneous  region  by  using  the  information  contained  in  scattered  electromagnetic  waves. 

Here  we  investigate  the  connection  between  two  relatively  diverse  areas  of  mathematical  physics, 
namely,  one-dimensional  inverse  scattering  theory  and  the  theory  of  almost  periodic  functions.  Our 
eventual  objective  is  to  determine  an  area  of  commonality  and,  more  specifically,  to  discover  inverse 
scattering  methods  which  can  be  applied  to  waves  scattered  or  reflected  by  almost  periodic  media.  This 
communication  serves  the  heuristic  purpose  of  investigating  the  feasibility  of  our  method,  first  sug¬ 
gested  in  (ll;  the  deeper  mathematical  questions  will  be  addressed  in  a  subsequent  paper. 

In  the  abstract  setting  of  Sabatier  [2),  inverse  scattering  theory  provides  the  mathematical 
methods  for  calculating  a  set  C  of  physical  parameters  for  a  set  S  of  scattering  data.  The  idealized  phy¬ 
sical  model  that  we  will  use  is  shown  in  Fig.  1.  In  this  idealized  model  the  set  S  consists  of  the  incident 
plane  wave  e*"  and  the  reflected  wave  r(Ar)e"'**,  where  k  is  the  free-space  wave  number.  The  set  C 
consists  of  /uq,  the  permeability  of  free  space,  and  €(k,  x),  the  permittivity.  The  inverse  scattering 
problem  here  is  to  determine  the  unknown  permittivity  e(A:.  jc)  from  the  reflection  coefficient  rik).  In 
general,  inverse  problems  start  with  limited  a  priori  knowledge  of  the  set  C  for  the  inhomogeneous 
region;  this  a  priori  knowledge  and  the  scattering  data  that  are  available  will  suggest  the  particular  inver¬ 
sion  method  that  should  be  used  and  the  physical  parameters  that  can  be  reconstructed.  In  contrast 
direct  problems  start  with  complete  knowledge  of  the  set  of  physical  parameters  and  seek  to  determine 
the  scattering  data.  Clearly  most  practical  problems  will  use  a  hybrid  approach.  For  example,  an 
inverse  problem  has  traditionally  been  solved  by  the  iterative  solution  of  the  relevant  direct  problem 
until  the  calculated  scattering  data  agree  with  the  experimental  data. 

The  mathematical  foundation  for  one-dimensional  inverse  scattering  theory  was  laid  by  the  Soviet 
mathematicians  Gel'fand,  Levitan  and  Marchenko,  circa  1 950  [3,4],  for  use  in  quantum  mechanical 
problems.  Subsequent  work  by  Kay  and  Moses  [5,7],  and  Newton  [6]  in  the  United  States  extended 
the  availability  of  these  results  and  expanded  the  theory  to  include  the  electromagnetic  problem.  Of 
particular  interest  is  the  treatment  by  Kay  (51  in  which  several  examples  display  the  usefulness  of  the 
Gel'fand-Levitan-Marchenko  theory.  Much  of  the  initial  work  [5.7,8.9,101  emphasized  exact,  closed- 
form  solutions  to  the  problem  when  the  scattered  data  were  represented  by  rational  functions.  Numeri¬ 
cal  [11.12,13]  and  approximate  [14,15]  solutions  provided  additional  understanding  of  this  problem.  A 
recent  summary  of  the  general  status  and  applications  of  several  inverse  scattering  methods  can  be 
found  in  a  recent  special  issue  of  the  IEEE  Transactions  on  Antennas  and  Propagation  [16]  and  in  a 
future  issue  of  Radio  Science  [17], 
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The  theory  of  almost  periodic  functions  has  its  origins  in  the  work  of  Bohr  [18]  and  Besicovitch 
[19],  performed  during  the  I920’s  and  I930’s.  Important  contributions  were  subsequently  made  by 
Wiener  [20]  and  others  [21],  More  recently,  these  functions  have  been  applied  to  direct  problems  in 
mechanics  [22,23]  and  electromagnetic  wave  propagation  [24,25,26,27].  These  functions  are  useful  in 
describing  certain  types  of  imperfections  in  periodic  structures  and  so  may  play  a  role  in  solid  state  phy¬ 
sics  and  in  the  design  of  optical  and  microwave  devices- 

Here  we  investigate  several  methods  to  determine  the  profile  of  a  structure  where  we  have  the 
a  priori  knowledge  that  this  structure  is  described  by  an  almost  periodic  function.  Trubowitz  [28]  has 
presented  the  inverse  scattering  theory  for  periodic  potentials.  After  reviewing  the  necessary  facts 
needed  from  inverse  scattering  theory,  the  theory  of  almost  periodic  functions  and  reflection  from 
almost  periodic  media,  we  use  several  examples  to  illustrate  the  application  of  inverse  scattering  theory 
to  reflection  from  almost  periodic  media.  Finally,  we  give  a  brief  discussion  and  several  general  conclu¬ 
sions  resulting  from  our  investigation. 

II.  MATHEMATICAL  BACKGROUND 


A.  Inverse  Scattering 


In  a  source-free  inhomogeneous  region  the  electric  and  magnetic  field  vectors  £  and  H  satisfy  the 
time- harmonic  differential  equations 


£  -I-  V 


V€ 


=  0, 


'2  //  +  Zl  X  V  X  //  -I-  e//  -  0, 


(1) 

(2) 


In  the  physical  model  of  Fig.  1,  we  have  assumed  that  the  plane-polarized  electromagnetic  field  is  nor¬ 
mally  incident  on  a  semi-infinite  inhomogeneous  region  whose  permittivity  f(k,  x)  is  assumed  to  obey 
the  dispersion  relation 


€(k.  x) 


(x  >  0) 


«o 


(JC  <  0) 


(3) 


where  €o 's  th®  permittivity  of  free  space,  so  that  the  behavior  of  the  time-harmonic  transverse  electric 
wave  amplitude  uik,  x)  is  described  by  the  differential  equation 


dx^ 


uik.  x)  + 


>  -  ir 


u(k,  x)  -  0. 


(4) 


The  transverse  magnetic  field  case  can  be  transformed  to  a  differential  equation  with  the  same  func¬ 
tional  form.  Here  q(x)  is  a  profile  function  related  to  the  density  of  the  unknown  region,  where  qix) 
is  assumed  to  be  piecewise  continuous.  Then  r(k)  will  be  an  analytic  function  of  k  so  that  the  conser¬ 
vation  of  energy  condition  is  expressed  as 

|r(/t)P  <  I.  (5) 


A  physical  description  of  the  inverse  scattering  procedure  [7]  is  obtained  by  using  the  Fourier 
transform  of  the  differential  Eq.  (4), 

U(x.  f)  -  Uix.  t)  -  q(x)  U(x,  l).  (6) 

9r 
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The  incident  plane  wave  is  represented  as  an  incident  impulse 


i/j-cU.  t)  =  S(x  -  t) 


which  produces  a  reflected  transient 


R(x  +  t)~  4-  r”  dk  -  i  T  (8) 

2ir  f*  ‘ 

where  the  in  the  second  term  are  the  residues  at  the  poles,  if  any,  of  r(k)  on  the  positive  imaginary 
/c-axis.  Because  of  causality,  we  must  have 

/?(jr  +  f)  -  0  ix  <  t)  (9) 

i.e.,  a  reflected  transient  is  not  produced  until  the  incident  pulse  has  interacted  with  the  medium.  It  is 
possible  (see,  e.g.,  Kay  (5,7])  to  relate  the  wave  amplitude,  U{x,  t),  in  the  inhomogeneous  region, 
X  >  0,  with  the  wave  amplitude  in  the  free  space  region,  U^ix,  t).  x  ^  0,  by  the  transformation 


where 


U(x,  t)  =  £/o(x,  t)  +  J  ^  Kix,  z)  Vq  (z,  t)dz 


Uq(x.  i)  =  8(x  —  r)  +  /J(x  +  f)  (x  <  0). 


From  physical  considerations  we  know  that  U(x,  r)  is  a  right-moving  transient  that  is. 

U(x.  f)  =  0  (x  >  /)  (11) 

so  that  U(x,  t)  depends  only  on  the  medium  traversed  by  Uoix,  t)  up  to  time  t.  Therefore, 
Kix,  z)  -  0  for  2  >  X.  Substituting  the  expression  for  the  field  in  the  free-space,  Uoix,  t),  gives  the 
integral  equation 

Rix,  t)  +  Kix,  t)  +  Kix,  z)  Riz  +  t)  dz  >^0  ix  >  t)  (12) 

J-y 

which  is  Kay's  version  of  the  Gel’fand-Levitan-Marchenko  integral  equation  and  is  to  be  solved  for  the 
function  Kix,  t).  Conditions  of  Kix,  t)  are  found  by  substituting  the  expression  (10)  for  Uix,  t)  in 
the  differential  (Eq.  (6));  the  function  Kix,  i)  satisfies  the  same  differential  equation  as  Uix,  f), 

Kix.  r)  -  '>  -  '>>  03) 

dx^  df^ 


subject  to  the  conditions 


Kix.  -x)  -  0 


—  Kix,  x)  •  2  qix).  (15) 

dx 

From  condition  (15)  we  see  that  the  inverse  scattering  problem  is  solved  for  qix)  if  we  can  solve  the 
integral  equation  (12)  that  relates  the  unknown  function  Kix.  t)  with  the  scattering  data  Rix  +  t). 

Thus  the  solution  of  the  one-dimensional  inverse  scattering  problem  that  we  are  investigating  has 
been  reduced  to  the  solution  of  Eq.  (12)  together  with  the  condition  (15).  Several  exact  and  approxi¬ 
mate  methods  of  solving  Eq.  (12)  have  been  reviewed  by  Jordan  [29];  the  most  important  for  obtaining 
closed-form  solutions  useful  for  engineering  applications  is  that  of  Kay  [7],  where  rik)  is  a  rational 
function  of  the  wavenumber  k. 


rik)  -  To 


ik  —  Hi)  ik  —  M:)  — 
ik  -  ki)  ik  -  ki)  — 
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where  k\,  *2.  ...  are  the  poles  and  mi.  M2.  •••  are  the  zeros  of  r(k)  in  the  complex  /r-plane.  Exam¬ 
ples  of  these  solutions  have  been  given  by  Jordan  et  al.  [9,10,11]  for  several  pole-zero  configurations; 
examples  of  numerical  and  approximate  solutions  have  been  given  by  Jaggard  et  al.  [13,15]  and 
Krueger  [12]. 


B.  Almost  Periodic  Functions 

In  the  restricted  sense  in  which  we  use  it  here,  an  almost  periodic  function  is  one  which  possesses 
a  discrete  or  line  spectrum.  This  implies  that  such  a  function  /(x)  can  be  written  in  a  generalized 
Fourier  series 

/(Jc)  -  Z  (17) 

n 

where  k„  is  the  wavenumber  or  spatial  frequency  of  the  mh  harmonic  and  a„  is  its  amplitude.  (Often 
we  will  refer  to  k„  as  the  tone  frequency  and  a„  as  the  tone  strength.)  This  equation  can  be  inverted  to 
produce  the  harmonic  amplitudes. 

o,  -  lim  C  f(x)e  dx.  (18) 

Clearly,  Eq.  (17)  represents  a  periodic  function  if  only  one  lone  is  present  or  if  all  of  the  tone  frequen¬ 
cies  are  commensurable.  In  this  case  Eq.  (17)  and  (18)  become  the  usual  expressions  for  the  Fourier 
series  and  Fourier  inversion  formula,  respectively. 

Although  almost  periodic  functions  do  not  possess  a  period,  as  noted  above,  they  do  possess  an 
almost  period  or  translation  number  r  such  that 

|/(x' -t- r)  - /(jt)|  <  rr,  .  (19) 

where  r  -  This  relation  provides  an  alternative  definition  to  the  more  restrictive  one  expressed 

by  Eq.  (17).  Another  property  of  interest  is  the  Parseval  relation  expressed  as 

lim  f  „  l/(x)|^  dx  -  ^  |a„|^  (20) 

These  and  other  characteristics  can  be  shown  from  the  defining  equations  (Eqs.  (17)  and  (19))  and  can 
be  found  in  the  literature  [18,19,20,23]. 

C.  Reflection  from  Almost  Periodic  Media 


If  the  potential  function  qix)  is  almost  periodic  an  approximate  theory  can  be  developed  for  the 
reflection  coefficient  using  the  method  developed  in  [26,27,30].  Here  we  examine  the  three  tone  case 
described  by 


q(x) 


Luai 
4 

0 


■y  COS  [(k  -I-  A)x]  -»-  COS  [kx]  -y  COS  [(k  -  A)x]  (Q  ^  v  ^  /) 


This  function  represents  a  modulated  periodic  structure  with  carrier  wavenumber  k  and  sideband 
wavenumbers  k  ±  1.  Clearly  A  is  the  tone  spacing.  Here  17  is  the  relative  amplitude  of  the  carrier  and 
m  is  the  modulation  index.  By  assuming  an  appropriate  form  for  the  electric  field  amplitude, 

M(x.  k)  -  fo(^)  e"*  BAx) 


Bo(x)  BAx)e  ^ 


(22) 
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(23) 


in  the  vicinity  of  the  primary  Bragg  resonance  defined  by 

k  =  k/2 


the  following  vector  coupled  mode  equations  (24)  result  under  suitable  approximations.  These  include 
small  or  moderate  amplitudes  (17  <  0.5,  ffi  <1)  and  relatively  closely  spaced  tones  (A  <  »c/4)  as 
discussed  in  [25,26]. 

f'  —  /■  8f  f  *  /  X  ^ 

-a' -  I  8b  J  -  I  x^f  (24) 


where 


and 


f  “ 


*  [Sq] 


B 


8b 


8-0  0 
0  80  0 

0  0  8. 


X  =  (x-  Xo  x^l 


80  «■  A/c  ""  Wc  —  k/2  8j.  *  80  ±  a 

X-  =  X*  *  Xo  “ 

The  prime  (')  indicates  differentiation  with  respect  to  x  and  the  dagger  (■*■)  denotes  the  Hermiiean 
transpose.  A  constant  loss  can  be  added  by  letting  80  and  become  complex. 


By  applying  the  appropriate  boundary  conditions  at  the  slab  ends. 


f  (0)  -  1 

[11 

£(/)  - 

Ifl 

/?_ 

0 

B(0)  - 

Bo 

-  R^*  BU)  - 

0 

B, 

0 

(25) 


the  problem  is  completely  specified.  Here  the  asterisk  (*)  denotes  complex  conjugate  and  the  last 
expression  in  Eq.  (25)  defines  the  reflection  matrix  R  while  the  second  expression  defines  the  scalar 
transmission  coefficient  f.  It  is  preferable  to  reformulate  the  above  two>point  boundary-value  problem, 
defined  by  Eqs.  (23)-(24),  into  an  initial-value  problem  for  the  reflection  matrix  to  avoid  problems  of 
numerical  instability  or  formidable  algebra.  This  leads  to  the  matrix  Riccati  equation 

R'  -  /’x  +  /8r  R  +  /R  8b  +  /  (R  x^  R  (26) 

with  initial  condition 


R(0)  -  10  0  01.  (27) 

This  formulation  is  particularly  amenable  to  numerical  solution  and  can  be  used  for  any  number  of 
tones.  A  result  for  the  three  tone  case  corresponding  to  (21)  is  shown  in  Fig.  2.  Here  the  scalar 
reflection  R  is  given  by 


R  -  RR* 


(28) 


(29) 


and  is  plotted  for  a  single  value  of  normalized  coupling 

xl  ”  +  Kx*  =  2 

as  a  function  of  normalized  frequency  /  =  21k/x  for  various  lone  spacings  and  losses.  Note  particu¬ 
larly  the  large  peak  at  the  primary  Bragg  resonance  defined  by  Eq.  (23),  or  /  =  0,  and  the  secondary 
Bragg  resonances  at 


it  =  (k  ±  A)/2. 


(30) 


These  results  suggest  the  following  observations: 

(1)  For  widely  spaced  tones  (bottom  plots  of  Fig.  2),  the  carrier  and  modulation  frequencies 
tend  to  act  independently.  Indeed,  this  can  be  shown  mathematically  since  the  vector 
coupled  mode  equations  separate  into  three  scalar  coupled  mode  equations. 

(2)  For  closely  spaced  tones  (top  plots  of  Fig.  2),  the  action  of  the  various  tones  is  no  longer 
independent  but  these  plots  are  suggestive  of  the  rational  function  reflection  coefficients 
displayed  in  |7,8,9]. 

For  these  reasons  it  is  natural  to  examine  these  two  classes  separately  in  the  inverse  problem.  For 
widely  spaced  tones,  use  of  resonance  theory  seems  to  be  the  most  promising  approach  while  for  closely 
spaced  tones  the  use  of  Gel’fand-Levitan-Marchenko  theory  appears  to  be  most  promising.  In  this  way 
the  t’.vo  methods  complement  each  other  so  that  a  wide  range  of  almost  periodic  media  may  be  charac¬ 
terized. 


III.  INVERSION  METHODS  FOR  ALMOST  PERIODIC  MEDIA 
A.  Widely  Spaced  Tones 

For  extremely  widely  spaced  tones  (A  »  \)  in  which  the  tones  act  independently,  the  results  of 
scalar  coupled  mode  theory  can  be  used  in  the  vicinity  of  each  Bragg  resonance.  For  example,  if  an  N 
tone  medium  is  defined  by 


q(x) 


1  ^ 
-■7  X 


(-1 


(0  <  X  ^  /) 

(x  <  0.  X  >  I) 


(31) 


where  «  KiOTKj  (;,  7-  I,  2,  ...  N),  then  the  scalar  reflection  R,  in  the  vicinity  of 

k  -  Kifl  is  given  by  131] 


IK,/ 


(32) 


/  coth  -v/xT"-^  /  -  i 

where 

V/  ”  Vt  (33) 

8,  /c  —  k,/2.  (34) 

The  peak  value  is 

l/Idmax  “  tanh  (v,/)  (35) 

which  can  be  inverted  to  get  ail  of  the  x,fs.  The  placement  of  the  peak  gives  the  values  of  the  k,'s. 
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If  the  combination  Xjl  is  not  desired,  but  instead  the  values  for  both  Xi  ^^d  /  are  to  be  deter¬ 
mined,  one  needs  to  examine  the  Bragg  resonances  at  the  sum  and  difference  wavenumbers  defined  by 

k  ^  -km  «|/2  ±  n  k-JT.  ±  p  kj/2  ±  ...  (36) 

where  m,  n,  p,  ...  are  positive  integers  or  zero.  As  the  simplest  example,  we  briefly  consider  the  two- 
tone  care  where  —  2  in  Eq.  (31).  Then  in  addition  to  the  values  of  K],  ki  and  xi  and  ^2  one  can 
look  at  the  sum  and  difference  resonances  which. occur  at  it  =  (k^  +  ki)I1  and  /c  =  1k2  —  k{J1  respec¬ 
tively.  The  peak  values  of  the  reflection  coefficient  here  are  determined  by  second  order  coupled  mode 


theory  as  described  in  [32].  The  result  is 

l^sumlmax  ”  tanh  iXsum^^ 

(37) 

and 

l^difflmax  “  tanh  (Xdifff^ 

(38) 

where 

Xsum  =  1^1  172  (*<1  +  '<2)/8 

(39) 

Xdiff  =  1?1  V2  -  »r2lV2  (ki  -I-  Ki). 

(40) 

By  using  these  values  of  xO'  Xsum/^ind  ^<(1^/ from  Eqs.  (32),  (37)  and  (38),  the  values  of  t}2. 
and  /  can  be  found  using  Eqs.  (33),  (39)  and  (40).  In  this  way  the  original  potential  function  q(x)  has 
been  reconstructed  up  to  a  relative  phase  between  periodicities. 


Referring  to  Tig.  3  the  steps  for  reconstructing  the  two  tone  potential  are  as  follows; 

(i)  Find  the  primary  maxima  evidenced  by  large  reflections  and  evaluate  ki.  «2,  )(i/and  xil 
from  the  reflection  coefficient  plot  and  Eq.  (35). 

(ii)  Find  the  secondary  maxima  at  /c  =  |k2  ±  kiI/2  and  find  Xsum/  and  xml  from  Eqs.  (37)- 
(38).  Calculate  tji,  tj2  and  /  from  Eqs.  (33),  (39)  and  (40). 

The  method  can  be  extended  to  the  three-tone  modulated  periodicity  of  (21)  and  to  potentials  with  an 
arbitrary  number  of  tones. 

The  criterion  for  the  applicability  of  this  method  is  that  the  tones  are  separated  sufficiently  so  that 
they  act  independently.  This  occurs  when  the  bandgaps  due  to  each  periodicity  are  spaced  in  frequency 
so  that  they  do  not  coalesce.  This  condition  is  described  by  the  relative  balance  of  the  periodicity 
amplitudes  with  the  wavenumbers  as  expressed  by 

|tj,|2  «  16  Ik,  -  K,|  (/,  J  ~  I,  2.  ...  N,  i  ^  j).  (41) 


B.  Closely  Spaced  Tones 

If  the  tones  are  closely  spaced  so  that  condition  (41)  does  not  hold,  then  the  simplified  expression 
(Eq.  (32))  for  the  reflection  coefficient  cannot  be  used.  The  traditional  approach  of  iterating  between 
the  direct  problem,  via  the  matrix  Riccati  equation,  and  the  experimental  data  leads  to  formidable  alge¬ 
bra  and  a  problem  of  physical  interpretation.  Alternatively,  the  inverse  problem  for  closely  spaced 
tones  can  be  approached  by  observing  that  the  graph  of  |r(/c)l^  for  a  rational  reflection  coefficient  with 
zeros  as  well  as  poles  resembles  the  behavior  of  the  scalar  reflection  for  reflection  from  an  almost- 
periodic  medium  with  three  closely  spaced  tones  (Fig.  2,  second  plot,  8  -  2x,  /•/  -  0.4).  The  band¬ 
pass  characteristics  of  the  almost-periodic  reflection  coefficient  in  the  vicinity  of  the  central  Bragg  reso¬ 
nance  suggests  that  poles  on  the  unit  circle  (Butterworth  poles)  could  be  used  as  a  first  approximation. 
The  number  of  "sidelobes"  of  l/-(/c)P  will  be  determined  by  the  number  of  zeros.  The  rapid  damping 
of  kl/clp  as  /c  —  oc  can  be  accounted  for  by  absorption  in  the  inhomogeneous  region. 


Combining  these  observations  suggests  that  a  rational  reflection  coefficient  with  the  form 

r(k)  =  _ 

Ml  ik  —  ki)  (k  —  ki)  {.k  —  kj) 


(42) 


where  the  locations  of  the  poles,  k\,  ki,  k},  and  the  zero.  mi.  in  the  complex  /t-plane  will  be  chosen  to 
meet  the  requirement  for  energy  conservation,  Eq.  (5),  as  well  as  approximating  the  almost-periodic 
reflection  coefficient.  This  particular  choice  for  rik)  was  obtained  by  combining  the  results  of  the  solu¬ 
tions  to  two  different  inverse  scattering  problems;  namely,  the  case  of  the  reflection  coefficient.  rs(Jt). 
due  to  a  complex  potential  function  (33)  and  the  case  of  a  2-pole  Butterworth  function  reflection 
coefficient,  rgik)  18],  Thus  we  write  r(k)  as 


where 


and 


rik)  =  tgik)  ■  r^ik) 
kik2 


(43) 


rgik) 


ik  -  ki)  ik  -  fcj) 


ky  /r-Mi 
rAk)  =  — 


Ml  k  -  ki' 


Here  we  set 


V  2  , ,  .. 

—  (I  -  /) 


la 


so  that 


k*  +  I 


Ml  *  Ha  -t-  S)  (0  <  S  «  a) 
k^  +  ia  +  5)^ 


(44) 


ia  +  8)^  k^ 

which  has  the  general  characteristics  of  the  plot  of  being  considered.  These  analytical  results  will  aid 
in  the  solution  of  the  inverse  problem  for  almost  periodic  media;  this  physical  model  can  be  interpreted 
as  an  inhomogeneous  medium  whose  effective  permittivity  has  dispersive  and  dissipative  characteristics 
which  we  now  determine. 

Applying  the  analysis  of  [33]  for  potentials  with  weak  absorption,  we  write 


'0 

- 

qix)  -  (fgix)  -t-  -1-  ikpAx)]  cos  kx 

(45) 

where  gix)  is  the  effective  complex  profile  function  for  almost  periodic  medium. 
[81  on  the  rgik)  of  (43),  we  find 

Using  the  method  of 

(I4-V2x)^ 

(46) 

using  the  method  of  [33]  on  r^k),  we  find 

<7s(x)  -  2  |/?,(2x)l'  [F(x)]  +  3  [/?,(2x)l2  [f(x)l^ 

(47) 

pAx)  -  RA2x)  lFi2x)V 

(48) 

where 

• 

a  8(1  -  8)  (x  ^  0) 

- 

ftAx)  - 

(49) 

0  (x  <  0) 
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(50) 


Fix)  =  tl  -  8  2 


The  inverse  scattering  method  has  been  demonstrated  with  a  reflection  coefficient  with  three  poles 
and  one  zero.  More  complicated  pole-zero  configuration  can  be  used  to  obtain  closer  agreement  with 
the  experimental  |/•(Ar)p  and  to  model  profiles  with  more  tones.  Losses  in  the  dielectric  region  are 
readily  included  in  the  inversion  method  by  poles  and  zeros  on  the  positive  imaginary  axis. 

IV.  DISCUSSION 

Inverse  scattering  theory  has  been  used  to  determine  the  permittivity  profile  of  a  one-dimensional 
inhomogeneous  medium  by  analyzing  the  properties  of  the  reflection  coefficient  amplitude  as  a  function 
of  incident  wavenumber.  The  problem  of  scattering  from  an  almost  periodic  medium  can  be  analyzed 
by  considering  two  cases  which  are  dependent  on  the  relative  spacing  of  the  tones. 

If  the  reflection  coefficient  is  composed  of  several  widely  spaced  primary  maxima  then  the  reso¬ 
nance  method  for  widely  spaced  tones  should  be  used.  The  number  and  spacing  of  the  tones  can  be 
determined  from  the  number  and  spacing  of  the  reflection  maxima.  A  check  on  these  tone  characteris¬ 
tics  can  be  made  by  examining  secondary  maxima.  The  length  and  amplitude  of  the  various  perturba¬ 
tions  can  be  found  from  the  heights  of  the  various  relative  maxima  of  the  reflection  data. 

If  the  reflection  coefficient  is  composed  of  closely  spaced  maxima,  then  the  complementary  inver¬ 
sion  method  must  be  used.  Here  the  reflection  coefficient  is  approximated  by  a  rational  function  of 
wavenumber  with  respect  to  the  principal  n^axima.  The  number  of  poles  and  zeros  of  this  reflection 
coefficient  function  and  their  location  in  complex  wavenumber  space  are  varied  until  this  function 
approximates  the  reflection  data.  Then  using  the  inversion  algorithm  of  Kay,  based  on  the  work  of 
GelTand,  Levitan  and  Marchenko,  the  potential  or  permittivity  is  reconstructed. 

Future  work  will  be  concerned  with  the  more  accurate  approximation  of  reflection  data  and  with 
finding  the  sensitivity  of  these  methods  to  reflection  data  errors. 
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IDEALIZED  PHYSICAL  MODEL  One  dimensional  scaiiering  of  a  lime-harmonic  wave 
of  unit  amplitude  from  a  half  space  characterized  by  tik,  .\  > 
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reflection  ^  as  a  function  of  normalized  frequency  ./  for  vari¬ 
ous  tone  spacings  S  for  a  three-tone  slab  of  length  /  as 
described  by  (2!)  The  tone  spacings  are  S  «  0.  2x-  b)(, 

and  8x  top  to  bottom  (sec  text!  and  the  relative  loss  Is  LI. 
Therefore,  the  top  two  plots  represent  closely  spaced  tones 
while  the  bottom  two  plots  display  tones  of  moderate  or  wide 
spacing.  The  modulation  is  constant  (m  —  I)  and  >(/ —  2 
Adapted  from  (27|. 
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SCALAR  REFLECTION  |R| 


Fig  3  —  TWO  TONE  REFLECTION  COEFFICIENT.  For  moderate  or  widely  spaced  tones  with 
wavenumber  k,  and  kj,  primary  Bragg  resonances  occurs  at  k  -  K|/2  and  k^/I.  Secondary  resonance 
occurs  at  subharmonic  and  superharmonic  wavenumbers  as  shown. 


